The B Meson Decay Constant from Unquenched Lattice QCD 
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We present determinations of the B meson decay constant /s and of the ratio fss/fB using 
the MILC collaboration unquenched gauge configurations which include three flavors of light sea 
quarks. The mass of one of the sea quarks is kept around the strange quark mass, and we explore 
a range in masses for the two lighter sea quarks down to ms/8. The heavy b quark is simulated 
using Nonrelativistic QCD, and both the valence and sea light quarks are represented by the highly 
improved (AsqTad) staggered quark action. The good chiral properties of the latter action allow 
for a much smoother chiral extrapolation to physical up and down quarks than has been possible in 
the past. We flnd fg = 216(9) (19) (4) (6) MeV and fsjfs = 1.20(3)(1). 

PACS numbers: 12.38.Gc, 13.20.Fc, 13.20.He 
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Accurate determination of the CKM matrix of the 
Standard Model and tests of its consistency and unitar- 
ity constitute an important part of current research in 
experimental and theoretical particle physics. Experi- 
mental studies of neutral Bd — Bd mixing, carried out 
as part of this program, are now well established and 
the mass difference AM^ is known with high precision 
Q. Uncertainty in our present knowledge of the CKM 
matrix element \Vtd\ is hence dominated by theoretical 
uncertainties, the most important of which are errors in 
Jb ^/Bb, where Jb is the B meson decay constant and Bb 
its bag parameter. Lattice QCD allows for first principles 
calculation of the hadronic matrix elements that lead to 
/b and /bVBb and in recent years the onus of reducing 
theoretical errors in determinations of \Vtd\ has been on 
the Lattice QCD community. In this article we address 
and significantly improve upon two of the errors that 
have plagued /s calculations on the lattice in the past, 
namely uncertainties due to lack of correct vacuum polar- 
ization in the simulations and errors due to chiral extrap- 
olations to physical up and down quarks. The generation 
of unquenched gauge configurations by the MILC collab- 
oration 2], which include effects of vacuum polarization 
from the strange plus two lighter dynamical quarks, has 
led to successful and realistic full QCD calculations of 
a variety of quantities involving both heavy and light 
quarks (a mi H, 111113 . Here we also take advantage 
of these well tested configurations. Another innovation 
in recent years has been to use the same improved stag- 
gered light quark action 10], which is being employed for 
sea light quarks and for light hadron physics, also for the 
valence light quarks inside heavy- light mesons This 
has been crucial for allowing heavy-light simulations close 
to the real world. Chiral extrapolation requirements are 
now much milder than in the past thus reducing effects 



coming from this source of uncertainty. 

In this study we work mainly with four of the "coarse" 
MILC ensembles with lattice spacing a around 0.12fm. 
These have dynamical light quark masses (in units of the 
strange quark mass) of mf/rris = 0.125,0.175,0.25 and 
0.5. We have also accumulated results on two of MILC's 
"fine" lattices with a ^ 0.087fm. On the fine lattices we 
use staggered valence light propagators created by the 
Fermilab collaboration. The heavy b quark is simulated 
using the same NRQCD action employed in recent studies 
of the T system '9j . For many of the ensembles the lattice 
spacing was determined from the T 25" — IS* splitting. For 
two ensembles where T results are not available, we used 
the heavy quark potential variable ri measured by the 
MILC collaboration 0,0. The bare s and b quark masses 
have been fixed by the Kaon and T masses, respectively 
0, and based on studies of light quark masses in 1^ we 
take as the physical chiral limit the point nis/niq — 21 A. 

The basic quantity that needs to be calculated in de- 
cay constant determinations is the matrix element of the 
heavy-light axial vector current between the B meson 
state and the hadronic vacuum. Taking, as is custom- 
ary, the temporal component of the axial current, in Eu- 
clidean space and in the B rest frame one has 
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In the last couple of years we have made considerable 
progress in reducing statistical errors in numerical deter- 
minations of this matrix element. We have developed 
better operators to create the B meson state on the lat- 
tice and fit to a matrix of correlators with different smear- 
ings. Details of smearings and matrix fits are similar to 
those in the T spectroscopy studies of reference and 
will not be repeated here. 

Table I summarizes results for the quantity <^q = 
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FIG. 1: The ratio ^<s> = ^s/^q versus mq/nis. The fuh hne 
through the data shows a fit to full QCD Staggered xPT (see 
text). Errors are statistical errors only. The fine lattice points 
were not included in the fit. The vertical line at niq/ms = 
1/27 A denotes the physical chiral limit. 



Jb^ ^JMb^ , where Bq denotes a " B" meson with a 
light valence quark of mass vriq. In the third column 

we show a?^^ , the result for <^q in lattice units 
when only the zeroth order lattice heavy-light current 
7570 is used. The next column shows 
q3/2<^^^ our results after one-loop matching and inclu- 
sion of \/M currents. All corrections to the heavy-light 
current at 0{Kqcd/M), 0(a,), 0{aa,), 0{aJ{aM)) 
and 0{as Aqcd/M) have been included. The dimension 
4 current corrections that enter into the matching at this 
order have been discussed in [l^ . The one-loop pertur- 
bative matching coefficients specific to the actions used 
in this study are given in 13]. One sees that the dif- 
ference between and $g is small, about 2 ^ 4% on 
the coarse lattices and ~ 7% on the fine lattices. The 
very small change on the coarse lattices may be partially 
accidental. There is cancellation between the 0{as) cor- 
rection to the zeroth order current and the 1/M correc- 
tions. The coefficient of the 0{as) term switches sign 
as one goes from a bare b quark mass of qMq = 2.8 on 
the coarse lattices to oMq = 1.95 on the fine lattices, so 
that the cancellation does not occur on the latter. In the 
last column of Table I we give results for in GeV^^^. 
The first errors are statistical and the second come from 
lattice spacing uncertainties. One sees that for most en- 
sembles scale uncertainties dominate over statistical er- 
rors. The scales, a~^, employed here are, in order of the 
most chiral to the least chiral ensembles, 1.623(32)GeV, 
1.622(32)GeV, 1.596(30)GeV and 1.605(29)GeV, respec- 
tively on the four coarse lattices and 2.258(32)GeV and 
2.312(31)GeV on the two fine lattices. 



TABLE I: Simulation results for $, = /b, a/A^s, ■ Sea 
(valence) quark masses are denoted by m/ (mq) and uq = 
[piag]^/* is the link variable used by the MILC collaboration 
in their normalisation of quark masses. See text for defini- 
tions of the last three columns. The second error in the last 
column comes from uncertainties in the scale a~^^^. 
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TABLE II: Simulation results for = <I>s/$g without and 
with 1/M plus one-loop corrections. 
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Table II shows results for the ratio £,<s, = This 
quantity, unlike $q itself, is not affected directly by errors 
in the lattice spacing. Several other systematic errors 
inherent in fg determinations, that will be discussed in 
more detail below, are also cancelled to a large extent in 
the ratio. For instance, one sees that going from ratios 
of to ratios of <I>'s that include 1/Ad and one-loop 
matching corrections, produces almost no change at all. 
The data for ^$ are plotted in Fig.l as a function of 
niq/ms. The full curve comes from a fit to formulas of 



staggered chiral perturbation theory (SxPT) 0,0,0 
and represents the prediction for full QCD. The vertical 
line at small ruq corresponds to the physical chiral limit 
niq/ms = 1/27.4. 

SxPT for heavy-light decay constants has been devel- 
oped by Aubin & Bernard in reference . For their 
formula reads, 



Co (1 + Aq 



analytic). 
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The term encompassing the chiral logarithms, = 
Sfsq/ (167r'^/^), is given in 0| and includes 0{a'^) effects 
coming from taste symmetry breaking, both in the mass 
splittings among light-light pseudoscalars and in lattice 
artifact hairpin diagrams. For the ratio ^$ we use the 
ansatz, 



= 1 + (A, 
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k 



amsf. (3) 



Nk was increased until the fit result for f<i, at 

mq/ms = 1/27 and its error had stabilized (in practice 
Nk = 2 was sufficient). Other ansaetze such as the direct 

ratio, IXtlXttZ'X^TlXZZl is the sea strange 

quark mass which, on the coarse lattices, is slightly larger 
than the true strange quark mass we use for valence 
strange quarks) or simple linear fits without any chi- 
ral logarithms were also tried as were fits with all the 
0{a^) taste breaking terms turned off. All these differ- 
ent chiral extrapolations lead to values for that 
differ at most by 3%. We fit simultaneously to the six 
coarse lattice points, 4 full QCD and 2 partially quenched 
(PQQCD) points, using fuU QCD and PQQCD SxPT 
formulas respectively. Fig.l shows just the full QCD 
curve. 

The terms A^ involve the BB*tt coupling gBn which 
is not known experimentally. We have carried out fits 
at several fixed values for g"^^ between g^^ — and 
g^^ = 0.75. Good fits were obtained {x^/dof w 1 or less) 
for gg^ < 0.5 with ^^^^''■^ differing again by less than 3% 



in the range ^, 
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— 1.21 ^ 1.24. We have also let gsrr 



float as one of the fit parameters and find g^^ = 0.0(2) to- 
gether with ^^P'^y'-^ = 1.21(2). This fit result for with 
the large uncertainty of Ag^^ — 0.2 shows that our data 
is not able to determine g^^ with any accuracy, the same 
message we get from the fixed gsTr fits, where a range of 
g%^ between zero and ^ 2 x Ag^^ all give acceptable 
fits. Fortunately, within this range ^2'''*'''''' is not very 
sensitive to gj^^. We take as our central value for ^ ^ 
the result from the floating gsTt fit, which we consider the 
least biased fit. This fit gives the curve shown on Fig.l. 
We then take ±0.03 as the error due to statistics and 
chiral extrapolation uncertainties, and which also covers 
the spread we observe upon trying different ansaetze and 
different ways of handling g^^. Remaining errors such 
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FIG. 2: $5 versus mq/rris. Errors include both statistical 
and scale uncertainty errors. The fine lattice points were not 
included in the fit. 



as those due to discretization and relativistic corrections 
and higher order operator matchings not yet included, 
will affect fs and fs^ in similar ways and largely cancel 
in the ratio. One expects their effects to come in at the 
level of the corresponding error in $q times a{ms — niq). 
We have already seen that 1/M and one- loop matching 
corrections cancel almost completely in Furthermore 
the two full QCD fine lattice points in Fig.l fall nicely 
on the full QCD SxPT curve fixed by the coarse lattice 
points indicating that any residual discretization errors 
in ^$ are smaller than the current statistical errors. Tak- 
ing all these arguments into account, we estimate a '--^ 1% 
further uncertainty in from these other sources. Our 

final result for fsj fs = 



is then 

Mbs 



/sJ/b = 1.20(3)(1) 



(4) 



We emphasize that the reason the chiral extrapolation 
errors are small here is because the light quark action 
employed in this study allowed us to go down as low as 
TOs/8 and only a modest extrapolation to the physical 
chiral limit was required. This differs from the case with 
Wilson type light quarks, where simulations have typi- 
cally been restricted to mq/ms > 0.5, i.e. to the region 
to the right of the heaviest data point in Fig.l. 

Fig. 2 shows the data points for $g itself for mq/ms < 
0.5 together with a full QCD SxPT fit curve. For chiral 
extrapolation of <I>g wc use directly eq.® with analytic 
terms ci(2m/ + mgd) + C2mq. We again carry out simul- 
taneous fits to the coarse lattice full QCD and PQQCD 
points. Fits with the coupling g^^ held fixed between 
0.0 and 0.6 afl lead to good fits with varying by 

4%. Allowing this coupling to float gives g'^^ = 0.1(5), 
which is consistent with the flxed gsir fit results, and 
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i^iphys.) ^ 0.496(20) GeV3/2 with again a 4% error. We 
take the 4% to be our best estimate for the combined 
error from statistics, chiral extrapolation and determina- 
tion of a-i. The full QCD SxPT curve in Fig.2 comes 
from the floating g^^ fit. We turn next to estimates of 
the other systematic errors in <J)(p''1'* ) . 

A major source of systematic error in i^'^P^v^-) is higher 
order matching of the heavy-light current. Although the 
one-loop contributions turned out to be small (as de- 
scribed above), in fact much smaller than a naive esti- 
mate of 0{aa) ~ 30%, we have no argument guarantee- 
ing this to be true at higher orders. Hence we allow for an 
©(a^) « 9% systematic matching error. This will be the 
dominant systematic error in our decay constant deter- 
mination. Another source of systematic error comes from 
discretization effects. The fine lattice points in Fig.2 lie 
about 3 ~ 5% lower than those from the coarse lattices. 
Since the statistical plus scale uncertainty errors on all 
our points range between 2 ~ 3%, it is not obvious how 
much of this difference comes from discretization effects. 
The size of fluctuations between independent coarse en- 
sembles is comparable to this difference. It should also 
be noted that the difference between the coarse and fine 
lattice data would disappear if it were not for the one- 
loop matching corrections (recall the 2 ^ 4% corrections 
on the coarse lattices versus the ^ 7% corrections on the 
fine lattices giving a 3 ^ 5% difference in the radiative 
corrections on the two lattices). In other words it is dif- 
ficult to disentangle discretization errors from radiative 
corrections. One could quote a combined discretization 
and higher order matching error again at the ^ 9% level. 
We opt instead to keep the 9% as the pure (and dominat- 
ing) ©(ttg) error and use a conventional naive estimate of 
©(a^tts) ~ 2% for discretization errors. As the last non- 
trivial systematic error we estimate uncertainties from 
relativistic corrections and tuning of the b quark mass |Q| 
to be at the ~ 3% level. Putting all this together we ob- 
tain ^(pi^y^-) = 0.496(20)(45)(10)(15) GcY^'^. This leads 
to our result for the B meson decay constant of 

/B = 0.216(9)(19)(4)(6)GeV. (5) 

The errors, from left to right, come from statistics plus 
scale plus chiral extrapolations, higher order matching, 
discretization, and relativistic corrections plus nib tun- 
ing respectively. Combining this result with our re- 
sult for /b./Zb, eq.O, one finds /b, = 0.259(32) GeV. 
This is very consistent with the direct calculation of 
fs^ published earlier in where we quote a value of 
0.260(29) GeV. 

To summarize, we have completed a determination of 
the B meson decay constant in full (unquenched) QCD. 
Our main results are given in eqns. and (j^l. The 
use of a highly improved light quark action has led to 



good control over the chiral extrapolation to physical up 
and down quarks. Better smearings have significantly 
reduced statistical errors. For the ratio JeJIb these im- 
provements translate into an accurate final result with 
errors at the ~ 3% level. For /b itself other systematic 
errors not yet addressed in the present study dominate 
and the current total error is at the ~ 10% level. The 
main remaining source of uncertainty comes from higher 
order operator matching. More studies should also be 
carried out on the fine lattices and on even finer lattices 
currently being created by the MILC collaboration, to re- 
duce discretization uncertainties. Errors in the scale 
need to come down for all the ensembles. Improvements 
on all these fronts are underway. Calculations of the bag 
parameter Bb have also been initiated. 
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